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Abstract 

We calculate analytically the asymptotic form of quasi-normal frequencies for massive scalar 
perturbations of large black holes in AdSs. We solve the wave equation, which reduces to a Heun 
equation, perturbatively and calculate the wave function as an expansion in 1/m, where m is the 
mass of the mode. The zeroth-order results agree with expressions derived by numerical analysis. 
We also calculate the first-order corrections to frequencies explicitly and show that they are of 
order 1/m. 
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Quasi-normal modes of black holes in asymptotically AdS space-times have received a lot of 
attention recently [1-17] and should shed some light on the AdS/CFT correspondence. They are 
derived as complex eigenvalues corresponding to a wave equation which is solved subject to the 
conditions that the flux be ingoing at the horizon and the wave-function vanish at the boundary 
of AdS space. The wave equation reduces to a Hypergeometric equation in three dimensions 
(AdSs) and can therefore be solved exactly [7, 13]. In higher dimensions, an analytic solution is 
not readily available, as the wave equation develops unphysical singularities. It has been possible 
to obtain numerical values for quasi-normal frequencies [2, 18-20]. 

More recently, we proposed an analytic method of calculating massless quasi-normal modes of 
large AdS black holes. We apphed the method to the calculation of low-lying frequencies [21] as 
well as in the asymptotic regime [22] in AdSs. The method was based on a perturbative expansion 
of the wave equation, which reduced to the Heun equation [23] , in the dimensionless parameter 
uj/Th, where uj is the frequency of the mode and Th is the (high) Hawking temperature of the 
black hole. This is a non-trivial expansion, for the dependence of the wave- function on uj/Th 
changes as one moves from the asymptotic boundary of AdS space to the horizon of the black 
hole. We showed that our results were in agreement with numerical results [2, 18]. 

Here we extend the discussion of [21,22] to the case of massive scalar modes. The quasi- 
normal frequencies determine the poles of the retarded correlation functions of dual operators in 
finite-temperature J\f — 4 SU {N) SYM theory in the large- A?^, large 't Hooft coupling limit [24]. 
They have also been recently studied in [25] where they were shown to arise in complexified 
geodesies. We shall derive a systematic analytic expansion of the frequencies in powers of 1/m, 
where m is the mass of the mode. Numerical results indicate [24] that in the limit of large mass, 

^r^2n{±l-i){n + h^-l) , n = l,2,... (1) 

where /i± = 1 ± Vl + ra^R'^- Our results confirm analytically the above asymptotic result of 
numerical analysis. Our method is an extension of the discussion in [22] where the numerical 
result (1) was confirmed analytically in the massless case (m = 0). We also calculate the first- 
order correction and show that it is of order l//i+. A similar procedure was followed in [21]. 
The metric of a five-dimensional AdS black hole may be written as 

ds^ = -/(r) de + ^^+ r^dnl , /(r) = ^ + 1 - ^ (2) 
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where R is the AdS radius and M is the mass of the black hole. For a large black hole, the metric 
simplifies to 



dr 



The Hawking temperature is 



where rh is the radius of the horizon, 



ttR^ 



1/4 



The wave equation for a massive scalar of mass m is 



-drir^ h{r) dr^) - 



(3) 



(4) 



(5) 



(6) 



r2/i(r) * 

We arc interested in solving this equation for a wave which is ingoing at the horizon and vanishes 
at infinity. These boundary conditions yield a discrete set of complex frequencies (quasi-normal 
modes). The solution may be written as 



(7) 



Upon changing the coordinate r to y. 



the wave equation becomes 



r 



4 y2_i 

where we have introduced the dimensionless variables 



- ^ - m^y 1 * = 



UJ 



m \p\ 



ttRTh' 



rh — mR 



(8) 



(9) 



(10) 



Two independent solutions are obtained by examining the behavior near the horizon {y ^ 1), 



*± - (y - 1)^'"/' 



(11) 



where is outgoing and is ingoing. We ought to choose for quasi-normal modes. 
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A different set of linearly independent solutions is obtained by studying the behavior at large 
r {y ^ oo). We obtain 



^ - y-''^ , h± = l± Vl + m2 (12) 

For quasi-normal modes, we are interested in the solution which vanishes at the boundary, hence 
it behaves as ~ y~'^+ as ^ oo. Combining this with the requirement that ^ — ^_ (eq. (11)) 
leads to a discrete spectrum of quasi-normal frequencies. 

By considering the other (unphysical) singularity at y = — 1, we obtain yet another set of 
linearly independent wave-functions behaving as 

^ - (y + 1)^*^/^ (13) 

near y = —1. There is no restriction on the behavior of the wave-function at this singularity. 

Nevertheless, as we shall see, it will be advantageous to isolate the behavior near y = —1. 
Isolating the behavior at the two singularities y = ±1, we shall write the wave- function as 

^iy) = {y-ir-/\y + l)-''/'F{y) (14) 

The selection of the exponent at the y — —1 singularity is arbitrary and not guided by a physical 
principle. We shall show that it leads to a convenient perturbative calculation of half of the 
modes. Selecting the other exponent, +oj/4, similarly yields the other set of modes without 
additional computational difficulties. The latter have the same imaginary part as the former, 
but opposite real parts, as will become evident. 

It is easily deduced from eqs. (9) and (14) that the function F(y) satisfies the Heun equation 



y{y' - l)i^"+ { (3 - + i^c^y - l} F' 



+ \2 (,T"^"V^"'^^ + ^^"'^4 " 4/ ^ ^ ° ^^^^ 
We wish to solve this equation in a region in the complex y-plane containing \y\ > 1, which 
includes the physical regime r > Vh- For large a), the constant terms in the respective Polynomial 
coefficients of F' and F in (15) are small compared with the other terms, so they may be dropped. 
Eq. (15) may then be approximated by the Hypergeometric equation 



in the asymptotic limit of large frequencies a). Two linearly independent solutions of (16) are 

J^i — F{a+,a-;c; —x) , — x^'''^F{l + a+ — c,l + a- — c;2 — c; —x) (17) 



where 

a±^h±-^u; , c^l-^co , (18) 

The solution of (16) which is regular at the horizon {x = 0) is the Hypergeometric function J-'i 
(eq. (17)). Matching the scaling behavior (12) at the boundary of AdSs, we obtain a different 
set of linearly independent solutions of (16), 

/C± = (x + l)-»±F(a±, c - a^; a± - + 1; l/(x + 1)) (19) 

We have /C± ~ x~"'^ as x ^ 00, and correspondingly, ^ ~ x~'^^ (from eq. (14) using (18)). 
Therefore, our choice is /C+, since it leads to ^ as x — > 00. It is a hnear combination of J^i 
and J^2 (eq. (17)). We easily obtain from standard Hypergeometric identities, 

/C+ = AoJ^i + BoJ^2 (20) 

where 

_ r(l-c)r(l-a, + a+) _ r(c-l)r(l + a+-a_) 

° r(l-a_)r(l-c + a+) ' ° r(a+)r(c-a_) ^ ^ 

For the correct behavior at the horizon, we demand 

-Bo = (22) 

which leads to two different conditions, 

c-a_ = -n + l , n = l,2,... (23) 

or 

a+^-n + 1 , n= 1,2,... (24) 
We obtain a set of quasi- normal frequencies by solving (23), 

tD„ = -2(l+i)(n + /i+-|) , n = l,2,... (25) 
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Another set is obtained from (24), 

cD„ = 2(1 + - 1) , n=l,2,... (26) 

To leading order, the two sets of modes have the same imaginary parts but opposite real parts 
{ujn ~ 2(±1 — i){n + h+)). This is in agreement with numerical results [18] (see eq. (1)). To 
go beyond the zeroth-order approximation, it is computationally advantageous to work with the 
first set of modes (25), which is a consequence of (23). The other set of modes, with positive real 
part, can be obtained without additional effort by choosing the exponent +a;/4 at the unphysical 
singularity {y — —1), instead (eq. (14)). 

The above approximation /C+ (eq. (19)) to the exact wave- function satisfying the Heun equa- 
tion (16) may be used as the basis for a systematic calculation of corrections to the asymptotic 
form (25) of quasi-normal frequencies. To this end, let us write the Heun equation (16) as 

{no + ni)F^o (27) 

where 

«» = (^-^>5J5 + |— --(3-— -j!'|;^-l(T-i-'j-'" P») 

and 

The zeroth-order equation, 

HoFo = (30) 
is the approximation we discussed above: the Hypergeometric eq. (16) whose solution is 

Fo = /C+ (31) 

(eq. (19)). By treating Tii as a perturbation, we may expand the wave- function, 

F = Fo + Fi + . . . (32) 

and solve eq. (27) perturbatively. Corrections to the quasi-normal frequencies (25) may then be 
obtained once an explicit expression for the correction to Fq has been calculated. 
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We shall restrict attention to the vanishing momentum case, p — 0. Extension to the more 
general case is straightforward. We have 

f.(x) = x;_m/°°^±^-x;,mj[°°^^ (33) 

where /C± are the two linearly independent solutions (19) of eq. (30) and W is their Wronskian, 

W = (a+ - a_)x-"(x + 1 (34) 

To study the behavior near the horizon (x ^ 0) , we may first analytically continue the parameters 
so that the integrals in (33) remain finite as x — > 0. Alternatively, we may introduce a cutoff ^ 
and notice that the part which is cut. 



SFi 



(X) = X:_(.) I - I !^ (35) 



remains finite as x — > 0. This is most easily seen by replacing the set of solutions /C± with 
and J^2 (eq. (17)). Therefore, the analytic continuation of the parameters does not affect the 
singularity at the horizon. Writing 



Fi{x) r^Ai + Bix^-" {x 0) (36) 



we therefore deduce 



where 



r(c- i)r(i + a± -g^) 

T{a±)T[c — Ozp) 



On account of (23), we have /?+ = 0, therefore. 



00 







Hence to first order, the quasi-normal frequencies are obtained as solutions of 

Bo + Bi = (40) 

where Bq is given by (21). 



7 



Let us first solve the n — 1 case in order to calculate the first-order correction to a)i. In this 
case a- — c (eq. (23)), and so 

X:+ = Fo = (a; + !)-«+ (41) 



Also, /3_ = and 



2(2^TT) ^° ^^^^ 



It follows that 



Using 

f oo 



/ dx \ ^ = S(l + A, 1/ - A)F(1, l + X;l + u; -1) (44) 



we deduce 



^ r(a_ — l)r(l + 0,1 — a_)a_ ^ , , , 

= 2rK)(a,-a_) {f(l.l + «-;l + a.;-l)-F(l,l + a-;«,;-l)} 



= S(a_-l,l + a+-a_) F(2, a_ + 1; a+ + 2; -1) 

2(a+ — a_)a+(l + a+j 

= S(a_-l,l + a+-a_) |^ + ...| (45) 
where we used a± — {i ± l)h+ + . . . and expanded 

\ —1 /i i\22/ \ —3 

« \ /I 1 \ a 2; / _ a 



ni,«;ft.)=^l--j.j +... (46) 

for large a and /9. 
Setting 

c - a_ = e (47) 

we obtain from eq. (21) 

Bq = eS(a_ - 1, 1 + a+ - a_) + . . . (48) 
Therefore, demanding Bq + Bi = (eq. (40)), we obtain the first-order correction 
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and the corresponding frequency 



Next, we discuss the case n — 2. We have 



/C+ = Fo{x) = {x + |l ^"^—ri^ + 

1^ a+ — a_ + 1 J 



1 /-oo 

«--!)(«- -2) 7o t:;^ 



Working as before, we deduce 

2(-- -,v— -/^u 

where 

7o = (a+ - a_ + 3)(a+ - a_ + 1) , 71 = -a+(3(a+ - a_ + 1) + 4) 
a+((3a+ + l)(a+ - a_ + 1) + 2a+) a?_(a+ + 1) 

72 = — ] , 73 = —7 

a+ — a_ + 1 a+ — a_ + 1 

After integrating, we obtain 

1 ^ 

^1 ^ TT7 ^ V7fc5(a_,a+ -a_ + A; - l)F(l,a_;a+ + A;- 1; 

2(a_ - l)(a_ -2) ^ 

This is to be compared with Bq (eq. (21)). Setting 

c — a_ = — 1 + e 

we have 



r(a_-2)r(l + at-a-) 
TM ' + 



We may write 



= "°'"TrLV'"°'' 5: 7t a_; + * - 1; -1) 

where 

, r(a_|_) r(a4. — a_ + A; — 1) 

" r(a+ + A; - 1) r(a+-a_ + l) 

Exphcitly, 

''''^ ^ ^ ~'~ a+ — a_ ~'~ (a+ — a_)(a+ — a_ — 1)) ' ^ l ~'~ 



72 = 3a+ + IH ^ , 73 = -«+ (61) 

a+ — a_ + 1 

Notice that all coefficients are o(/i+), whereas we expect to finally get an expression which is 
o(l//i+). This means that terms have to conspire so that both o(/i+) and o(l) contributions 
cancel. In comparison, recall that in the n — 1 case, only o(l) terms appeared and they were 
seen to cancel. 

Let us expand the coefficients to the desired order, 

7^ = a+ + ^ + ^ + ... , 7 =-3a+ ^ 62 

72 = 3a+ H + . . . , 73 = -a+ (63) 

We may also use (46) to expand the Hypergeometric functions (higher-order terms are not needed, 
because J^^k — •^(1/^+))- We obtain 

, , a+ a-(ka+ + (k — 2)a_) , 

F 1, a_; a+ + k-l; -1 = + \^ ^ ^3 + " " " ^4 

a+ + a_ (a_|_ + a-)"* 

It is now a straightforward exercise to show that all terms o{h^) and o(l) cancel and we finally 

obtain 

_ r(a_ - 2)r(l + a+ - a_) f + 5a^a_ + 4a+a?. + 2a^ 



2r(a+) i (a2 - a2_)2 

r(a_-2)r(l + a+-a_) / 12i - 1 



2r(a+) t 16/iH 

Using eqs. (40), (57) and (65), eq. (56) reads 



+ ... (65) 



12i - 1 



The corresponding frequency is 



= + l + + (67) 

For a general n (eq. (23)), /C+ = Fq contains a Polynomial of order n — 1 in (x + 1)~^ (the 
Hypergeometric function in eq. (19)). Thus, the complexity of the calculation increases with 
increasing n. One obtains coefficients which are of order Even though they can be seen to 

cancel order- by-order, leaving an expression of order l/h+ after the dust settles, yielding 

cD„ = -2(l+i)(n + /i+-|) + o(l//i+) , n = l,2,... (68) 
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the calculation of the integral (39) is uninspiring. We calculated the coefficient of in (68) 
exphcitly in the cases n = 1 (eq. (50)) and n = 2 (eq. (67)). It would be desirable to have an 
expression for the first-order correction (39) which is manifestly of order 

It would also be interesting to extend the above results to more general black holes in AdS^ and 
compare with existing numerical results (see, e.g., [17]). The number of (unphysical) singularities 
increases if d 7^ 3, 5, so the analysis is expected to be more complicated. 
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